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Perturbation-Minimizing Frequency Assignment in a 
Changing TDMA/FDMA Cellular Environment 

Dong-Wan Tcha. June-Hyuk Kwon, Taek-Jin Choi, and Se-Hyun Oh 


Abstract —A robust operation of assigning frequencies to 
requirements in a time-division multiple-access (TDMA) [or 
frequency-division multiple-access (FDMA)] cellular system 
should cope with environmental changes such as short-term 
demand rises and long-term capacity expansions while always 
keeping the required realignment process as simple as possible. 
In this paper, we consider the so-called perturbation-minimizing 
frequency assignment problem (PMFAP), the objective of which is 
to assign available frequencies for newly generated requirements 
with the minimum change in the existing frequency assignments 
while meeting the interference-related constraints. For PMFAP, 
we propose a heuristic algorithm based on the so-called Bv-Dw 
rearrangement technique, which can also be applied to the classic 
frequency assignment problem (FAP) with a slight modification. 
Two kinds of computational experiments, one for each of the 
above two problems, are performed to demonstrate the powerful 
features of the proposed solution method not only in its suitability 
for real-world frequency management, but also in solving the FAP. 

Index Terms —Assignment with rearrangement, frequency as¬ 
signment problem, TDMA/FDMA. 


I. Introduction 

T HE RECENT demand explosion for mobile commu¬ 
nication services, together with the finite frequency 
spectrum allocated to this service, makes the issue of frequency 
assignment ever more important for the design and operation 
of such systems. The cellular system is founded on frequency 
reuse between cells sufficiently apart from each other, so that 
it must prohibit use of the same and/or adjacent frequency 
channels (or simply frequencies) in each cell and its neigh¬ 
boring cells to prevent radio interference. Three kinds of such 
interference-related constraints—cochannel, adjacent channel, 
and cosite—have to be satisfied for the cellular system (see [1] 
and [2]), which are aggregately called here the compatibility 
constraints. The objective of frequency assignment problem 
(FAP) is to assign frequencies for given channel requirements 
with minimum frequency span while meeting the compatibility 
constraints. Though the FAP is originally for frequency-divi¬ 
sion multiple-access (FDMA) systems, it is also well defined 
and very critical for a time-division multiple-access (TDMA) 
system with multiple carriers (see [3]). 
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However, the FAP is not quite suited for being directly used 
by the system operator which has to meet real fluctuations 
and temporal increases in demand with a given number of 
available frequencies. Consider the case when the final solution 
of the FAP is so different from the preassigned pattern of the 
operative system as to require a full-scale frequency reshuffle. 
The reshuffling would heavily burden the system both in time 
and in cost. From the system operator’s point of view, the 
smaller the change in the preassignment, the better, particularly 
for a short-term fluctuation in the traffic environment. Another 
point hampering its direct use arises from doubting that avail¬ 
able frequencies beyond the frequency span obtained from 
solving the FAP are of no value. This realization motivates us 
to consider the so-called perturbation-minimizing frequency 
assignment problem (PMFAP), the objective of which is to 
meet newly generated channel requirements with a given 
number of available frequencies while keeping the number of 
frequency reassignments as low as possible. 

The problem of assigning frequencies for cells with newly 
generated requirements was studied by Lee et al. [4]. They con¬ 
sidered, however, the model with only cochannel interference 
constraints. Furthermore, their model does not allow frequen¬ 
cies assigned to the existing requirements to be altered, which 
certainly reduces the flexibility required to accommodate more 
channel requirements. 

In Section II, we briefly review the classic model of FAP, 
along with two representative solution strategies. Taking the 
notational framework, we then present our PMFAP formula¬ 
tion. In Section III, a heuristic algorithm is proposed, each it¬ 
eration of which consists of three operations, requirement se¬ 
lection (RS), unforced assignment (UA), and forced assignment 
with rearrangement (FAR). The proposed algorithm for PMFAP 
is shown to effectively solve the FAP itself by always attempting 
to assign to a requirement under consideration in the algorithmic 
process the lowest numbered frequency in the given set of or¬ 
dered available frequencies. The version for FAP is summarized 
in Section IV. In Section V, we report two kinds of computa¬ 
tional experiments. The first experiment is to show the effective¬ 
ness of our algorithm for PMFAP in minimizing the number of 
frequency reassignments. The focus of the second one is placed 
on demonstrating the competitiveness of the version of our algo¬ 
rithm modified for FAP in solving the FAP itself. Finally, Sec¬ 
tion VI provides concluding remarks. 

II. Notation and Problem Statement 

Frequencies are assumed to be evenly spaced, so they can be 
identified with positive integers 1, 2, 3, etc. We denote the uni¬ 
versal set of requirements by X — {aq, X 2 , • • •, x n }, to each 
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member of which a frequency is assigned such that the frequen¬ 
cies assigned to two requirements, for example, x,. X j £ X , 
should be separated by at least r,j. Traditionally, the whole re¬ 
quirements are grouped into cells, based on which the compat¬ 
ibility constraints are defined with indexes referring to cells. 
However, for exposition brevity, we let the compatibility con¬ 
straints be represented by compatibility matrix C = {c,y} as 
defined above with indexes referring to requirements. Then in¬ 
teger Cjj denotes the cosite constraint when two requirements 
Xj and xj belong to the same cell, whereas both adjacent and 
cochannel constraints correspond to the case where two require¬ 
ments belong to different cells. 

Given set X and n X n matrix C, a frequency assignment is 
represented by n vector 1 — (fj), where component /, repre¬ 
sents the frequency assigned to requirement x, £ X. For fu¬ 
ture use, a two-tuple (x,. fj ) will be called a requirement-fre¬ 
quency pair (or a r-f pair) and the set {(x,. fj)\xj £ X} a re¬ 
quirement-frequency match (RFM). A frequency assignment F 
is called a feasible (frequency) assignment, or an RFM is fea¬ 
sible, if it satisfies the following condition: 

I fi ~ fj | > Cjj, for all Xj, Xj £ X and i f j. 

Then, the objective of FAP is to find a feasible assignment F 
such that 


Then, the objective of PMFAP is to find a feasible frequency 
assignment F = (fj) or a feasible RFM for X such that 

P(F):= Y, I&) 

xi ex° 

(i.e., the level of perturbation on the initial assignment) 
is as small as possible. 

Throughout this study, we shall employ the requirement-fre¬ 
quency matching rule to assign the lowest numbered frequency 
whenever possible. Under this rule, it would be instructive to 
distinguish between two cases, one with a sufficiently large 
number of available frequencies and the other not. In the former 
case, an algorithm for PMFAP would proceed to minimize 
the associated frequency span, which also makes itself fit 
for solving an FAP. Consider the latter case where there are 
not enough available frequencies to meet the compatibility 
constraints of the requirement set X, as often encountered in 
practice. The natural goal of this type of PMFAP would be to 
minimize the number of requirements which remain unmatched 
and then to reduce the level of the above perturbation. This 
preemptive preference of two objectives is well answered by 
the algorithmic process that all the available frequencies are 
always orderly and exhaustively searched to obtain a feasible 
r-f pair for each requirement. 


S(F ) := max fj (i.e., the required frequency span) 
is as small as possible. 

The FAP is NP-complete, since it can be viewed as a gener¬ 
alization of the NP-complete graph coloring problem [1], [6]. 
The traditional solution approach for FAP has thus been on de¬ 
veloping heuristics by extending the graph coloring algorithms 
[2], [5], [7]—[9]. Most of such FAP heuristics are based on se¬ 
quential assignment of frequencies according to some heuristic 
order. That is, the typical strategy is to rank whole requirements 
in descending order of assignment difficulty and then to as¬ 
sign frequencies one by one to the ordered requirements via a 
number of differing assignment rules. The paper of Sivarajan et 
al. [5] may be the representative of such studies. They defined 
the degree of a cell and used it as a measure on the difficulty 
of assigning frequencies to requirements in the cell. They fur¬ 
ther proposed four requirement ordering methods based upon 
this measure and used two representative assignment strategies 
which are frequency exhaustive strategy (FES) and requirement 
exhaustive strategy (RES). 

For PMFAP, we denote the set of available frequencies as 
L — {1, 2, • • •, (}. LetX 0 withm = |X°| be the set of initially 
given requirements for each component Xj of which /° has been 
assigned. Given such an RFM 0 = {(%i, /, 0 )|x, £ X 0 }, we aim 
at finding a feasible assignment Fora feasible RFM for a set 
X expanded from X°. 

For a feasible assignment F of X, we define, for every x, £ 

X°, 


I(xj ) 


1, if fit ft 

0, otherwise. 


III. Heuristic Algorithm for PMFAP 

In this section, we propose a heuristic algorithm for PMFAP. 
Each iteration of our algorithm consists of three operations: 
requirement selection (RS), unforced assignment (UA), and 
forced assignment with rearrangement (FAR). Let Q (C X) 
be the set of requirements to each member of which a cer¬ 
tain frequency is already assigned during the algorithmic 
process, including the preassigned requirement set X°, and let 
Q = X\Q. A requirement-frequency pair, (x-j, fj) for Xj £ Q, 
is said to be feasible to the present RFM if it does not conflict 
with the RFM in terms of the compatibility constraints, i.e., 

I fi ~fj\> Cjj, for all (xj, fj) £ RFM. 

For the purpose of explaining the three operations, it will be 
convenient to introduce the following definition. 

Definition 1: Consider a subset of the RFM, the constituent 
r-f pairs of which conflict with the r-f pair (xj, fj), Xj £ Q. 
We shall call the subset of requirements in Q, corresponding 
to the pair, as the fj blocker of Xj, denoted as B(xj, fj). So 
if the frequency assignments for requirement set B (x, ■ fj) are 
undone, then fj becomes a feasible frequency for or the r-f 
pair (xj, fj) becomes feasible, with the remaining pairs of the 
RFM. 

If B(xj, fj) — 0 for fj £ L, then fj is immediately as¬ 
signable to Xj, symbolized as fj — x,. The case of nonempty 
B(xj, f) is denoted by fj^Xj. 

The RS operation selects a requirement which is most diffi¬ 
cult to assign a frequency among the requirements in Q, and 
this difficulty measure is a generalization of the saturation de¬ 
gree proposed by Brelaz [10] for GCP. For this, the assignment 
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difficulty of an unassigned requirement x, (6 ()) is defined as 
the number of infeasible frequencies with respect to the RFM, 
that is, the number of /, ’ s such that /, — . 

Once an unassigned requirement x, is selected, the UA op¬ 
eration finds the lowest frequency in the list L feasible to the 
RFM. That is, if an (6 L) such that B(xj , /, ) = 0 is found, a 
successful assignment is made to proceed to the next iteration. 
Otherwise, the FAR operation takes over. 

It is the FAR operation that makes the proposed algorithm 
prominent. Assume that for an Xj £ Q , there is no with -a 
.r, with respect to the present RFM. Then, the FAR attempts to 
assign a frequency in L to requirement x; with the minimum 
enforcing overhead, i.e., with the minimum perturbation on the 
present RFM. The essence of FAR is to identify a subset Rlx,) 
of Q, which is defined as follows. 

Definition 2: The R(.r,) for Xi (6 Q) is defined as a minimal 
subset of Q, each requirement of which can be simultaneously 
reassigned with an alternative feasible frequency, so that the x, 
can be assigned a frequency feasible to the realigned RFM. 

There can be a number of R{xfy s for an Xj (£ Q) from its 
definition as a minimal subset of (). To identify one such R(xf), 
we examine a sequence of /, blockers for /, 6 L until a termi¬ 
nation criterion is met. Each time a B (./■,. /,) of Xj. is generated, 
we undo the corresponding portion of frequency assignment in 
the RFM and try to assign an alternate feasible frequency to 
each requirement of B (x, • /,) by the UA operation. If these 
frequency reassignments are successfully made, B'-.ig.. /,) be¬ 
comes R(xi) by itself. 

In case such a frequency reassignment cannot be made for 
some requirement xj in B(xj, f,) , one proceeds to identify 
B (Xj . f j) and attempts to reassign an alternate feasible fre¬ 
quency to each xu £ B(xj, fj). Such B(xj , fj)’s are blockers 
at the second depth level or simply 2-DL. Generalizing this, 
we define the kth depth-level (fc-DL) blockers for the r-f pair 
(-G• fi). - ! 'i £ Q- 

Now suppose that fc-DL blockers have been successively 
identified up to level l to make way for assigning /, to Xj 
(G Q). Let Bixp. f p ) be a nonempty /-DL blocker. The task 
is then to check whether each requirement of the blocker 
is reassignable with an alternate feasible frequency. For the 
systematic checking of (l + 1)-DL blockers, we define the 
requirement-frequency fixing, RFF((7t:,', /,:), [x p , f p )), to be 
the set of r-f pairs, which have been successively matched 
along the vertical path from (x,. f,) down to (x p , f p ) in the 
RFF search tree. Thus, when at the (7 + 1)-DL, we identify 
those blockers of r-f pairs which do not conflict with those 
pairs in RFF((a;,-, /, ). (x p . f p )). Note that the cardinality of 
RFFff.r,. /,), (.Xp. f p )) is not less than its path length, as seen 
in the following example. 

Now with the notation defined above, the FAR operation is 
illustrated via the following simple example in Fig. 1. 

In this example, the RFM is {(a, 3), (6, 2), (c, 1), 
( d , 2), (e, 3), (/, 1)} and the list of available frequencies L is 
{1, 2, 3}. The existence of cochannel interference between a 
pair of requirements is represented by an edge connecting the 
corresponding two nodes on the figure. 

We shall first elaborate on the vertical (depth-first) search 
of our frequency assignment process for requirement <j. 


3 2 



To assign a feasible frequency to < 7 , we first generate 
three blockers: B(g, 1) = {c}, B(g, 2 ) = {b, d}, and 
B(g, 3) = {e}. Focusing on requirement c in the first 
blocker above, we move on to the 2-DL with the updated 
RLM:= {(a, 3), (6, 2), ( d , 2), (e, 3), (/, 1), (g, 1)}. Since 
(c, 1) conflicts with RLL((#, 1), (g, 1)) := {(g, 1)}, 

we consider only two blockers: B(c , 2) and B{c , 3). 
Lrom B(c , 2) = {of} and 3 — 7, the associated ver¬ 
tical search is fathomed, giving rise to an Rig) — {c, d}. 
Lrom B(c , 3) = {a, e} and 2 -a e, we move on to 

the 3-DL for an unassigned requirement, a, by setting 
RLL((< 7 , 1), (c, 3)) := {((/, 1), (c, 3), (e, 2)}. Note that 

the inclusion of (e, 2) in the RLE is to simplify the search, 
which in effect implies the nonexhaustiveness of our search. 
Lor requirement a, we have two pairs, (a, 1), and (a, 2), not 
conflicting with the RLE. Generation of 3-DL blockers shows 
B(a , 1) = {/} and B(a, 2) = { 6 }. Searches along this 
vertical direction continue in this manner. 

Going back to the original requirement g, recall that we have 
B(g, 2) = {6, d} and B(g, 3) = {e}. The UA operation shows 
1 —^ 6, 3 —^ c/, and 2 —1 e. Lor these two r-f pairs, ( <j. 2) and 
(fj. 3), two sets of i?(< 7 )’s are found as {6, d} and {e}, respec¬ 
tively. Note that we have another set of R(g), {c, d}, obtained 
as above, shown in Lig. 2. 

Lrom the above simple example, we have witnessed how 
often the process of checking the reassignability of an /, 
blocker, is repeatedly executed in generating an R(g). The 
following lemma shows the complexity of this process. 

Lemma 1: The problem of checking the reassignability of a 
nonempty blocker B (./■,. /,) is NP-complete. 

Proof: Suppose that the frequency assignments for 
B (x,. ff) are undone, rendering the REM reduced that much. 
Lor each requirement Xj £ B(xi, /,:), let L{xj) (C L) be the 
set of frequencies feasible to the reduced REM. Without loss of 
generality, let the compatibility constraint for any two require¬ 
ments in B{.r,. f,) be either zero or one. The checking problem 
under this restriction becomes a GCP on graph G — (V, E), 
where node vj corresponds to requirement xj in B(xj, ff), 
edge (vj, Vk) is defined if Cjk = 1, and colors of L(xj) are 
available for node Vj. The GCP is known to be NP-complete, 
and so is the checking problem. □ 

The PAR operation is encoded to render the so-called vth 
breadth-level and ruth depth-level (Bv-Dw) procedures. In the 
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Bv-Dw procedure, we consider blockers only within the car¬ 
dinality of v (i.e., | B(xi, fi )| < v) and limit the number of 
successive downward searches to w. This can be visualized via 
the associated RFF search tree: each node can have at most v 
dangling nodes and the path lengths are at most w. 

The reassignabilities of blockers at the same depth level are 
systematically checked in the increasing order of blocker car¬ 
dinality. The procedure starts from B 1, moves on to B2, and 
continues in this manner until Bv at the D 1 level and then re¬ 
peats this process for each unassigned requirement in blockers 
of the next lower level. However, the complexity of FAR opera¬ 
tion actually prohibits the direct implementation of this general 
Bv-Dw procedure. Viable alternative search strategies would 
be accentuating either of two tracks: breadth and depth. Two 
procedures, Br-D 1 and Bl-Dr, then arise. 

Procedure: Br-Dl procedure of the FAR operation. 

Step 0 ) Start with the given RFM° and an unassigned re¬ 
quirement Xi 6 (). 

Step 1) Find {B(xi, fi) | V fi 6 L} and set k 1. 

Step 2) If k > r, terminate (/* fail to produce a feasible 
assignment '/). Otherwise, let 

L k := {fi e L\\B( Xi , fi)\ = k}. 

Step 3) If L k — 0, set k k + 1 and go to Step 2). 
Otherwise, select a frequency f* in L k . 

Step 4) If all requirements in B (./■,. /*) are found re- 
assignable, update the RFM as found and stop. 
Otherwise, set L k L k \{f*} and go to Step 3). 

Considering the NP-completeness of the reassignability 
check, it is certainly not worthwhile to continue the search 
beyond a certain breadth level, say Br, though there is a slim 
chance of the reassignability of the remaining large blockers. 
At this point, it is instrumental to note that the search pattern 
not only saves the associated computational burden, but also 
conforms with our objective of minimizing the level of pertur¬ 
bation on the present RFM. 


For the other procedure, Bl-Dr, identifying blockers and 
checking their reassignability get more complex as the depth 
level increases. So instead of listing the Bl-Dr procedure, 
we elaborate on the B1-D2 procedure which limits the search 
within the D2 level. 

Procedure: B1-D2 procedure of the FAR operation. 

Step 0) Start with the given RFM 0 and an unassigned re¬ 
quirement Xi 6 (). 

Step 1 ) (D1 level) Call the Bl-Dl procedure. If x-i is as¬ 
signed, update RFM and stop. Otherwise, let 

L 1 := {fi e L\\B(x t , fi )| = 1}. 

Step 2) If L 1 = 0, terminate (/* fail to produce a feasible 
assignment */). Otherwise, select the frequency f* 
in L 1 and set RFF:= {(a:,-, /*)}. 

Step 3) (D2 level) For xj 6 B(xi, /*), let 

{/, C /. /T.r,,./M 1 

S:B(Xj. fj) not conflicting with RFF}. 

Step 4) If L 2 — 0, set L 1 := L 1 \{/*} and RFF:= 0, and 
go to Step 2). Otherwise, select the frequency f in 
Lf. 

Step 5) If the requirement in B(xj , f) is found re- 
assignable, update the RFM as found and stop. 
Otherwise, set L 2 L 2 \{f} and go to Step 4). 

IV. New Algorithm for FAP 

The FAR operation can be directly employed to solve the 
classic FAP, in which RFM 0 := 0 according to our setting. 
The proposed algorithm will be called the frequency exhaus¬ 
tive strategy with rearrangement (FESR), which are obtained by 
modifying the FES [2], [5], 

Recall that the FES for FAP proceeds as follows: starting at 
the top of the ordered requirement list, it consecutively assigns 
the lowest possible frequency which does not violate the com¬ 
patibility constraints with the assignments already made. Let S k 
denote the current frequency span after completing the assign¬ 
ment operations of step k in the FES process. At every step k, 
if the lowest feasible frequency f* is less than or equal to the 
current span S k -i, the FES can then assign a frequency without 
increasing the current span; otherwise, the increase of the span 
is inevitable in the FES. The former operation is called the fre¬ 
quency-insert and the latter the frequency-append. 

This FES is a straightforward implementation of the sequen¬ 
tial graph coloring algorithm [12]. The weakness of such an al¬ 
gorithmic process is the rigidity that r-f pairs, once matched, 
cannot be reshuffled later on. The FESR modifies the rigid FES 
by substituting the FAR operation for each frequency-append. A 
simple, but effective FESR algorithm is now presented, where 
the employed FAR operation is the Bl-Dl procedure. Listing 
of a general FESR algorithm incorporating a Bv-Dw procedure 
is also avoided for exposition brevity. 

Algorithm: Frequency exhaustive strategy with rearrange¬ 
ment. 

Step 0) Start with the ordered requirement list. Let k := 1, 
S k := 1, and RFM := 0. 



394 


IEEE TRANSACTIONS ON VEHICULAR TECHNOLOGY, VOL. 49, NO. 2, MARCH 2000 



Fig. 3. The 21-cell system. 

Step 1) Set k k + 1. If k < n, find the lowest frequency 
f* for Xu feasible to the RFM. Otherwise, stop (/* 
a solution is at hand */)• 

Step 2) (Frequency-Insert) If f* < update RFM by 

adding r-f pair . /*). Set .S';. := ,S';._i and go to 
Step 1). 

Step 3) ( Frequency-Rearrangement ) 

a) Let / := 1. 

b) If / < /*, obtain a set B(xu, /)■ Otherwise, go 
to Step 4). 

c) If \B(xu, f) | = 1, find the lowest frequency 
/ for x £ B(xu, f) that (x, f) is feasible to 
the reduced RFM and /). Otherwise, set 
f f + 1 and go to Step b). 

d) If f < f% update RFM by adding both (xu, f) 
and (x, /). Set Su max{Sfc_i, /, /} and 
go to Step 1). Otherwise, f f + 1 and go to 
Step b). 

Step 4) ( Frequency-Append) Update RFM by adding 

{xk, /*) and S k f*, and go to Step 1). 

Note that our term for the frequency-insert is the UA op¬ 
eration in PMFAR while our FAR is associated with the fre¬ 
quency-rearrangement. 

V. Computational Results 

The 21-cell system in [5] and [11], the configuration of 
which is given in Fig. 3, is still used to generate test problems. 
To demonstrate the powerful features of the proposed algo¬ 
rithm, we perform two kinds of computational experiments. 
The first experiment aims at showing the effectiveness of our 
algorithm for PMFAP in minimizing the level of perturbation, 
i.e., of requirement-frequency rematches for RFM 0 . Its effec¬ 
tiveness, however, may not be sufficiently illustrated due to the 
nonexistence of a study in the literature directly dealing with 
the PMFAP. Motivated by this realization, we perform another 
experiment, the focus of which is placed on demonstrating the 
superiority of FESR, the version of our algorithm modified for 
FAP, over the existing FAP algorithms in solving the FAP itself. 

A. Experiment on PMFAP 

We consider two scenarios on channel requirements as shown 
in Table I. The first one typifies the situation of short-term de¬ 
mand rises, whereas the second one a system expansion with 
newly added cells. The first column of each scenario shows the 
number of existing requirements to which frequencies are al¬ 
ready assigned, and the second one the increased number of ag¬ 
gregate requirements. 

As for the interference constraints in this experiment, we use, 
to facilitate the comparisons with the existing studies, the tra- 


TABLE I 

Channel Requirements for the 
21-Cell System 


Number of Requirements 


Cell 

Scenario 1 

Scenario 2 

existing 

new 

existing 

new 

1 

10 

12 

11 

14 

2 

12 

15 

9 

13 

3 

10 

14 

13 

13 

4 

11 

11 

10 

12 

5 

12 

13 

0 

9 

6 

12 

14 

0 

7 

7 

13 

16 

8 

11 

8 

11 

12 

9 

10 

9 

12 

17 

12 

14 

10 

10 

11 

10 

11 

11 

13 

15 

11 

13 

12 

14 

15 

0 

8 

13 

13 

16 

0 

9 

14 

15 

19 

14 

15 

15 

12 

13 

8 

10 

16 

12 

12 

12 

13 

17 

13 

17 

10 

11 

18 

10 

15 

11 

12 

19 

10 

12 

0 

6 

20 

15 

18 

13 

13 

21 

15 

17 

0 

8 

total 

255 

304 

161 

232 


ditional cell-based values for the adjacent channel and cosite 
constraints which are, respectively, set at two and five for sim¬ 
plicity. For the cochannel constraint characterized by N c , how¬ 
ever, we consider two values of N c , 12 and 7, to check the effect 
of the severity of interference environment. Note the differences 
of parameter values of constraints from those of our require¬ 
ment-based formulation of PMFAP. For example, i\> = 12 here 
implies that the same frequency can be reused between two cells 
with radius distance of six. The number of available frequencies 
l is also varied between two values to check the effect of fre¬ 
quency availability. 

With this input data specification. Tables II and III summa¬ 
rize the performance of our algorithm. Throughout this experi¬ 
ment, we set the initial requirement-frequency match, denoted 
by RFM 0 , as the best solution with the minimum span among 
those generated by the FAP algorithms listed in [5]. To demon¬ 
strate the superiority of solutions generated by our algorithm, 
two measures are used: the number of requirements which re¬ 
main unmatched even after running our algorithm, denoted by 
U(X\X°), and the number of matches altered from the initial 
RFM 0 , denoted by _P(RFM°). Table II shows an interesting phe¬ 
nomenon: the quality of solutions generated by the algorithm 
with B2-D1 is not noticeably better in both measures than the 
most basic algorithm with Bl-Dl, whereas the B1-D2 outper¬ 
forms the Bl-Dl though with extra computational burden. The 
plain performance of B2-D1 after the application of Bl-Dl 
may be reasoned as follows: the RFM obtained after the Bl-Dl 
rearrangement is so densely ordered that reshuffling the existing 
matches with a large cardinality may sometimes do harm on the 
existing compact match pattern. 

Table III is provided to additionally accentuate how well our 
algorithm serves its objective of minimizing the level of per¬ 
turbation. For the new aggregate demand of each scenario, we 
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TABLE II 

Comparisons Between Various Versions of Our Algorithm 





Algorithm with Bl-Dl 

Algorithm with B2-D1 

Algorithm with B1-D2 

Scenario 

Nf 


U(X \ X°) 3 > 

P(RFM°f'> 

CPU fl > 

U(X \ x°) 

P(RFM°) 

CPU 

U{X\X°) 

P(RFM °) 

CPU 


7 

130 

19 

7 

2.5 

19 

11 

3.1 

16 

20 

40.1 

i 


150 

0 

2 

0.3 

0 

2 

0.4 

0 

2 

0.3 


12 

160 

22 

11 

4.7 

23 

14 

5.1 

15 

42 

60.6 



180 

2 

4 

1.1 

2 

4 

1.1 

0 

12 

5.3 


7 

100 

14 

5 

1.4 

15 

7 

1.5 

12 

11 

10.8 

2 


120 

1 

1 

0.3 

0 

3 

0.3 

0 

3 

0.6 


12 

120 

32 

5 

3.3 

32 

5 

3.6 

32 

13 

82.2 



140 

4 

4 

0.8 

4 

4 

0.9 

9 

_2_ 

5 

18.7 


b cluster size 

2) number of available frequencies 

3 ) number of requirements which remain unmatched 
11 number of matches altered from RFM° 

b computation time in seconds 


TABLE HI 

Comparisons Between the Existing FAP and Our Algorithms 


Scenario 

N c 

Sivarajan et al. 

Algorithm with B2-D1 

Algorithm with B1-D2 

span* 

P(RFM°) 

U(X \ X°) 

P(RFM a ) 

U(X \ X°) 

P(RFM°) 

i 

7 

150 

227 

0 

2 

0 

2 


12 

181 

227 

1 

4 

0 

6 

2 

7 

124 

139 

0 

1 

0 

1 

_ 

12 

144 

135 

0 

_I_ 

0 

7 


* the minimum among the spans obtained by eight algorithms of Sivarajan et al. 


TABLE IV 

Performance Comparisons on FAP 


Problem 

LB 

Sivarajan et al. 

[5] 

Wang et al. 
[15] 

Sung et al. 
[16] 

FESR with 
Bl-Dl 

FESR with 
B1-D2 

i 

533 

533 

533 

533 

533 

533 

2 

533 

533 

533 

533 

533 

533 

3 

381 

381 

381 

381 

381 

381 

4 

427 

445 

433 

436 

434 

433 

5 

529 

529 

529 

529 

529 

529 

6 

309 

310 

309 

309 

309 

309 

7 

258 

270 

263 

268 

265 

260 


applied the well-known algorithm for FAP by Sivarajan et al. [5] 
to obtain the RFM. The first column shows the associated span, 
and the second the relatively large value of P(RFM°) indicating 
the full-scale reshuffling on the existing match. Taking this span 
as our l , the results of the algorithm with B2-D1 and B1-D2 
are listed in Table III. It is remarkable to find the outstanding 
records of two algorithms in both performance measures. 

B. Experiment on FAP 

Seven FAP examples widely used in the literature are again 
adopted to compare the performance of our FESR on solving 
the FAP with those of the well-known FAP algorithms [5], [15], 
[16], A version of the FESR is made for each of four require¬ 
ment ordering methods in [5]. Each performance record listed 
in Table IV is the best one among those obtained by the above 
four versions of the corresponding algorithm. The column of 
LB shows the tightest lower bound among those in [11], [13], 
and [14], 


The figures of resulting frequency spans in Table IV show that 
our FESR, particularly that with B1-D2, outperforms all the ex¬ 
isting FAP algorithms, even the most recent one with significant 
margin in some cases. This is indeed an impressive performance 
record of the FESR on solving the FAP when considering that 
the FESR is simply a modification of the algorithm for solving 
the PMFAP. 


VI. Concluding Remarks 

We addressed the so-called PMFAP, the objective of which 
is to meet newly generated channel requirements with a given 
number of available frequencies while keeping the level of per¬ 
turbation as low as possible. The PMFAP was defined to suit 
the real-world frequency assignment practices in TDMA and 
FDMA cellular systems where the level of perturbation is pre¬ 
ferred to be kept minimal for any demand fluctuation. 
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For PMFAP, we have developed a heuristic algorithm based 
on the Bv-Dw type of frequency rearrangement technique. The 
proposed algorithm, though originally developed for PMFAP, 
can be applied to the FAP with a slight modification. 

To demonstrate the powerful features of the proposed solu¬ 
tion methods, we perform two kinds of computational experi¬ 
ments. The first experiment is to show the effectiveness of our 
algorithm on solving the PMFAP, and the second one the com¬ 
petitiveness of the modified version on solving the classic FAP. 
This computational experience convinces us that the algorithms 
based on the Bv-Dw search not only are efficient enough to be 
directly used for real-world frequency management by system 
operators, but also outperform all the existing FAP algorithms 
even on solving the classic FAP itself. 
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